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Abstract 



04 ' We present a general scheme to realize the POVMs for the unambigu- 

^ , ous discrimination of quantum states. For any set of pure states it enables 

■ us to set up a feasible linear optical circuit to perform their optimal dis- 
crimination, if they are prepared as single-photon states. An example 
of unknown states discrimination is discussed as the illustration of the 
general scheme. 

^ ■ PACS number(s): 03.65.Bz, 03.67.-a, 03.65. Ta 

Discrimination of the members of a set of n quantum states {pi\ [i = 
1, 2, • • • , n) is a standard task in quantum communication protocols. Typically 
these states are nonorthogonal pure states or even mixed states, so it is impossi- 
(— I ■ ble to achieve one hundred percent correctness in the discrimination process. We 

\ can eliminate the chance of error but the price we will pay is that we must allow a 

^ ■ chance of inconclusiveness. This senario is called unambiguous state discrimina- 

tion (USD) T, and it has found applications in quantum crytographic protocol 
13 and quantum algorithm j3j, etc. The physical methods proposed to do USD 
k> \ include linear quantum optics 0j , ion trap architecture [^j and nuclear magnetic 

■ resonance [H]. Some experiments on the USD of nonorthogonal photon states 
' [7] have been realized thus far. To perform the USD of nonorthogonal states, 

we need to use general positive operator-valued measures (POVMs) instead of 
orthogonal projectors, and their realization in the original signal Hilbert space 
of the measured system is usually impossible. Through Neumark's theorem 
however, a POVM can be realized in the extended Hilbert spaces by perform- 
ing unitary transformations and von Neumann projections together, and then 
the physical implementation of this POVM is feasible due to the fact that any 
discrete finite dimensional unitary operator can be constructed in laboratory 
using optical devices [HI- In finding the necessary unitary operators, we used to 
require the success probabihties, pi = Tr{piIlk)Si^k, of the POVM elements Ilfc 
and the corresponding inconclusive probabilities, qi — 1 — pi jlOj . However, for 

* Department of Physics and Astronomy, Hunter College of the City University of New 
York, 695 Park Avenue, New York, NY 10021, contact email; bhe98@earthlink.net 



1 



the USD of a set of unknown states, we don't have the information available to 
obtain these probabilities and, therefore, are unable to find the required unitary 
transformations by the existing methods. 

In this letter, we present a general scheme to realize any POVM for USD 
only with its elements 11^ we have set up to unambiguously measure a particular 
set of pure or mixed inputs. The required unitary (orthogonal) transformations 
are found in an extended 27V dimensional Hilbert space if the dimension of 
the original signal Hilbert space is N. If the input states for discrimination 
are prepared as one-photon states, we can build a linear-optics circuit with at 
most lossless beam splitters to implement the POVM that unambiguously 
discriminates them. To reach the optimal performance of our set-up, we just 
need to choose the proper parameters (transmission and reflection coefficients) 
of the beam splitters such that the unitary (orthogonal) transformations in the 
extended space will achieve the same effect as the optimal POVM does in the 
original space. It would also be equally possible to construct this device for 
the quantum states in the forms of any other type of systems, such as elec- 
tron, neutron, atoms, etc. Moreover, this scheme is generalizable to the optical 
implementation of the Kraus operators of a POVM, and of any possible finite 
linear map given as a quantum operation, the general aspects of the physical 
realization for which has been discussed in 

The construction of the unitary transformations in the extended space by 
this general scheme is independent of the knowledge about the input states, so 
it can be applied to realize the USD of the unknown states, for which the inputs 
with certain symmetry are prepared with the copies of the unknown states as 
in |12l Take the inputs and the POVM in ^2] for example, two unknown 
qubits I'l/'i) and \ip2), which are randomly distributed unit vectors on the Bloch 
sphere, 

IV,) = cos(0,/2)|O) +sin(0,/2)e''^'|l) , (1) 
are prepared with their copies to produce the following quantum registers, 

I*,) - Hl)A\i^2)B\i'^)c, (2) 

where i — 1,2, with the probabilities of rji and 772 repectively, and the POVM 
elements to unambiguously discriminate them are constructed as follows: 

Ili = A\Ai = ki{P^'c®lA) 

Us = 4^2 = k2{PTc®lB), (3) 

where P]^^ and are the projectors into the anti-symmetric spaces of the 
respective digit locations. The coefficients fci and k2 can be properly tuned that 
this POVM will reach the optimum performance in the USD of these inputs. 
Meanwhile, the probabilities of inconclusiveness are determined by the operator, 

no = A5Ao = /-ni-n2. (4) 

Since the inputs are unknown to us, the only quantity to indicate how well this 
POVM performs is the average succss probability. For the above example of 
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1 copy of IV'i) and \tlJ2) used as the reference in the registers, the maximum 
average success probality is 1/6 if the preparation probabilities 771 and 772 are 
equal [T^ . 

Generally, by a direct sum, we append an M dimensional ancilla space A 
to the N dimensional original Hilbert space H to realize a POVM through 
Neumark's theorem. With quantum optics techniques (see M. Mohseni et 
al. in [7]), n arbitrary initial input states = Sjgjdijaj|0), where J = 

{1, 2, • • • , TV + M} and dtj = ioi N + 1 < j < N + M , are prepared as lin- 
ear combinations of single-photon states through multi-rail representation |14| . 
and are mapped by a unitary transformation in the extended space to the final 
states, the parts of which in the original Hilbert space H are orthogonal and can 
be distinguished among themselves by counting the photons in different output 
ports. The output photons recorded in the ancilla ports then correspond to the 
inconclusive results. 

We start our direct sum realization of all such POVMs with the observa- 
tion that the inconclusive operator XIq is positive and Hermitian in any of an 
orthonormal basis {|ei)}, where i = 1,2, ■ ■ ■ , N, in our N dimensional Hilbert 
space H. Then we will find a unitary transformation U that transforms the 
general orthonormal basis {le^)} to a unique orthonormal basis {la,)} (up to 
some permutation), where Hp is diagonalized: 



N 



C/HoiJt = ^Q|a,)(a,|. (5) 



1=1 



Because Hq is a positive operator and (0|Ho|0) < 1 for any |0) in H, all its 
eigenvalues satisfy < q < 1. From this fact we obtain the following well- 
defined operators: 



N 



Ao = Al = unY,V^M){'^^\]u, (6) 



and 



(7) 



if we represent them with the general orthonormal basis {|ei)}. With these 
operators we construct S and other three unitary, or more exactly orthogonal 
transformation, operators in the extended 2N dimensional space as follows: 



{I~AlAo)i -Ao 
Ao (I-AlAo^. 



(8) 



The other three such operators are obtained by putting the minus sign in the 
upper right sub-matrix of E to the other three sub-matrix blocks, respectively. 
It is straightforward to prove E^^S = SE^ = / with the operators defined in Eq. 
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(6) and Eq. (7). In the most general situation when we have n inputs to be 
unambiguously discriminated among themselves, I — AqAq in the square root 
equals Sjg/IIj, where 7 = {1, 2, • • • , n}. 

We take E to act on a set of states {pi} {i = l,2,---,n), which are to 
be distinguished between each other in a USD process, in the extended 2N 
dimensional Hilbert space: 

H = - 4^o)ip.(/ - AlAo)i {I - AlAo)ipiAl \ . . 

o J \ Aopiil - aIAo)-2 AoPiAl J' ^ ' 

where the blank blocks and the o sub-matrix represent the parts with all the 
entries 0. The trace of the upper left diagonal block gives the success probability 
Pi of unambiguously determining pi because Tr{pillj) = piSij for the Ilj's, and 
the trace of the lower right diagonal block gives the failure probability Qi in 
the ancilla space A, and in the whole extended space K = H (B A we have 
Pi + qi = Trp, = 1. 

If {pi} is a set of linearly independent pure states we will prove 

that the parts of their outputs after the action of S are mutually orthogonal 
in the original signal Hilbert space H. Before the transformation, they are 
extended to the inputs = Hi^i) , 0)'^ , where (0) represents a N-tuple of 

zero's, (0, 0, • • • , 0), in A. The output states are obtained as follows: 



{I-AlAo)^ -Ao \/|Vi> 
Ao {I-AlAo)i J\ 

{I - AlAo)-^\^i) 
Ao\^i) 



(10) 



Then the inner product of the outputs for any pair of different \ipi) and \tjjj) in 
H is 

n 

(V,-|(/-4^o)*(/-4^o)i|V'i) = {^j\Y,AlAkm = 0, (11) 

k=l 

and in A is 

(V,|4^o|V'i) = (12) 

For the mixed states we act E on the products piPj {i ^ j) in the same way 
as in Eq. (9), and find that the trace of the outputs in H vanishes while in A is 
the same as those of the inputs piPj. Therefore, the unitary (orthogonal) trans- 
formation S in the extended 2N dimensional Hilbert space realizes a scheme to 
unambiguously discriminate any set of quantum states {pi}. 

Next we need to find the way of how to optically realize the unitary transfor- 
mation S, given that the inputs are linear combinations of single-photon states. 
Because the representation of E may be complicated with the general orthonor- 
mal basis, we apply the following unitary transformation in K to reduce it to 
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be diagonalized in 4 blocks: 

U ^\( {I 



AlA 



Ao 



-Ao 
iI-AlAo)i 



I 



N 

J2 VI - Ci\at){ai\ 

N 

i=i 



N 
i=l 



N 



(13) 



We use Tp q to define an identity matrix with the entries Ipp, Ipq, Iqp and Iqq 
replaced by the corresponding 0(2) matrix elements. If we multiply the matrix 
obtained in Eq. (13) by Tp q, only the entries at {p,p), ip,q), {q,p) and {q,q) 
positions will be changed with all other entries intact. Taking out the entries 
on the ith and (i + iV)th rows and columns for example, and multiplying them 
on the right by the corresponding 0(2) sub- matrix in Ti_|_jv,i, we have 



VT^ - 
^ VT 



cos Ui 

sin 9j 



— sm ti 
cos 9i 



cos 9 i \/T^^Ci ~ sin 9 i^/ci — sin^^Vl 



sin^iVl — Ci + cos9iy/ci cos 9i\/l — Ci 



- COS( 

sin 9, 



Ci 

c~ 



(14) 



If we choose tan6'i = — y T^T" ^^'^ off-diagonal elements vanish, the diago- 

nal elements will be 1 and an identity matrix will be obtained. This two dimen- 
sional rotational transformation matrix, with its diagonal and off-diagonal ele- 
ments given as the transmission and reflection coefficient respectively, can be im- 
plemented by a lossless beam splitter, if the inputs are single-photon states. Sim- 
ilarly, if we right-multiply the whole matrix successively by T2n,n,T2N-i,n-i, 
■ ■ •, Tjv+i^i, we will obtain 



/ N 

VI - Ci\ai){ai 



N 



N 



\ 



J2 V^\ai){ai 



- J2 V^\ai){ai\ 

i=l 

N 

J2 VI - Ci\ai){ai 
1=1 



2N,NJ-2N-1,N-1 ■ • ■ r^v+ia 



/, (15) 



an identity matrix in the extended 2N dimensional space K . Therefore S rep- 
resented by the basis is given as 



E = 



Ar+1,1^ JV+2,2 



. • • • 



(16) 



and can be exactly implemented by N lossless beam splitters j^. Following the 
notation of P. Kok et al. in we choose the relative phase shift f = for 
these beam splitters. 

In most realistic applications, we also need a post-process unitary transfor- 
mation V after the action of E so that the output states, which have already 
been orthogonal, will be mapped to different final output ports for measure- 
ment. Together with the former procedure, the whole unitary transformation 
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Ut (in the extended space K) required to realize the USD of a set of inputs is 
constructed as follows: 



Ut = 



V 



V 



I N 

%/l - Ci\ai){ai\ 

1=1 

N 



N 

E Vc7|a»)(aj| 



N 



V 



i=l 1=1 



U 



U 



(17) 



which consists of the pre-process transformation U, S in the form of 4 di- 
agonalized sub-matrices and the post-process transformation V. We don't 
actually need U and V in the ancilla space A, so the total number of the 
beam splitters required to perform the USD of the inputs has an upper bound 
2 X (A^(A^ — 1) /2)+N = N^, a value determined by the dimension of the Hilbert 
space of the measured system. 

As an application of this general scheme, we show how to realize the USD 
of a pair of unknown qubits discussed in |12| . We assume that the preparation 
probabilities rji and 772 of the inputs are equal for simplicity (the unequal rji 's 
situation is realized in the same way but the matrix elements involved would 
look more complicated) and give the optical realization procedure to achieve 
their optimal USD. 

The quantum registers in Eq. (2) are expanded as — S^-dj where 
{\x)} is {|000), |001), • • • , |111)} and can be realized by 8 modes of single-photon 
state = a[.|0) (fc = 1,2, •••,8) respectively. With these basis vectors we 
perform the following transformations: 



IXii) = 



no) 
no) 



2' ' V 2' 

where z = 0, 1, to map them to a set of another orthonormal basis 
{\g,)} = {|000), |100), |r,oi), Ixoi), |011), hii), Ixii), Im)}- 



m) 



zOl), 



(18) 



(19) 



These orthogonal transformations can be implemented by beam splitters in the 
input states preparation period, and the advantage of doing this is that the 
multiplicity of the eigenvalues of IIq will be conveniently demonstrated if we 
represent it with this basis. Then the quantum registers for the USD process 
are represented with this basis as follows: 



I* 



(2)\ 



(20) 



where 



I* 



(1) 



: cos ■ 



■ sm 



cos ■ ^ 



■ sm 



1 

-^cos^^sin^e"^! 



cos 



cos 



92 



COS 



%/2cos^cos^sin^e*'^2 



^ cos -f 

cos^ ^ sin ^e' 

- -l=cos2 ^sin^e*'^^ 



(21) 
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(2) (2) 

represented by {|100), |?7oi), |xoi)}, and I^I/^ ') and l^'g ), which arc spanned 
by the set of basis vectors {|011), Ixii)}) just carry the above |*i^''),|*2^'') 
components but with one cos ^ in the terms replaced by sin ^e'"^* because there 



is one more digit 1 in the basis vectors. The ') and ') (resp. ) and 
|^'2^^)) in the joint space of and |^'2) are given as 



(0)\ 



cos 



sm 



|cos||000) 

— e'^vi sin— e'"^^ I 111), 
2 2 ' ' 



(22) 



and |^I'2''^), 1^*2'''') are obtained by interchanging 61, (f>i with 62, (p2 for all the 
factors in Eq. (22). We should bear in mind that we know nothing about 6iS 
and ipiS because and [^2) are randomly distributed on the Bloch sphere. 

With the basis {\gi)} the inconclusive operator IIo for the optimal USD 
measurement is given as the following direct sum: 



,(3)\ 



Eo = To e Ti e Ai e Aq, 



(23) 



where Fq and Aq are 1x1 matrices 1 acting on the spaces spanned by |000) 
and I 111), respectively, and 




(24) 



acting on the space spanned by {|100), I7701), lxoi)}j and Ai acting on the space 
spanned by {|011), |?7ii), |xii)} looks similar to Ji with the signs of some off- 
diagonal elements changed. There is a rotation. 




(25) 



which can be decomposed into three 2 dimensional rotations and is thus realized 
by three beam splitters together, to get Fi diagonalized: 



000 
I 
001 



(26) 



and Ai is diagonalized to the same matrix or has the same eigenvalues as Fi. 
Now we have the diagonalized operators. 



(/-F;)i = (7-A;)i = 




(27) 



7 



acting on the respective sub-spaces of H, and the unitary (orthogonal) operator 



/ 



V 



-1 



,(28) 



; 



acting on the extended space K. 

We therefore just need to look at the USD of \^[^^) and |*2^^) for the whole 
problem. The output parts of these states in the original Hilbert space after U 
and E are 



I* 



I* 



= (/-r;)^c/W|*('^^ = 
= (/-rl)^c/«|vi/«) = 



--i= cos |l sin ^e*"^! cos % + ^ cos^ ^ sin ^e*'^^ 
cos |l sin ^e'vi cos ^ + ^ cos^ ^ sin ^e"^^ 


^ cos ^ cos ^ sin ^e"f^ 
- ^ cos 1^ cos ^ sin ^e*'^^ 



cos 



2 6(2 



^sin^e^'^i 
__i.sin"fe^'^i cos'^ 



2 

,2 62 



(29) 







respectively. At this step they have been orthogonal to each other, and then we 
add a Hadamard gate, 



1 1 

%/2 ^2 



V2 





V2 









1 



(30) 



as the post-process to rotate these two states such that they will have non-zero 
components in different output ports. Thus the final states will be unambigu- 
ously distinguished between each other once the photons are recorded in different 
output ports. 



The process on l^-f^) and is similar too. I^i"''): l*2''0 and 1^-^' 

|\E'2 ) will be zero vectors after we map them with the corresponding transfor- 
mations to the final states. Together with 2 zero components produced in other 
sub-spaces like those in Eq. (29), there is a 4 dimensional joint space of {"^i) 
and 1^2): the vector components in which only contribute to the inconclusive 
result. The average of the inner products of the outpout vectors in Eq. (29) 

is 1/12. Putting it together with the corresponding contribution from |^]^ ') or 
(2) 

1^*2 ), we have the total average success probability 1/6 for the optimal USD 
of l^-i) and 1*2) [EI- 



,(2) 



(0)\ 



,(0)x 
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All the transformations in the process only need 14 beam splitters (including 
some totally reflecting mirrors) to implement. Moreover, in the whole process, 
we don't need to acquire any information about 9iS and ipi^s of the randomly 
distributed qubits. By the general scheme we propose in this paper, we realize 
the USD of the unknown states prepared as in Eq. (2). 

B. He would like to thank Dr. Z. Wang for helpful discussions. 
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